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Macroscopic Superconducting Current through a Sihcon Surface Reconstruction with 

Indium Adatoms: Si(lll)-(y7 x V3)-In 
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Macroscopic and robust supercurrents are observed by direct electron transport measurements 
on a silicon surface reconstruction with In adatoms (Si(lll)-(A/7 x y/3)-ln). The superconducting 
transition manifests itself as an emergence of the zero resistance state below 2.8 K. I — V char- 
acteristics exhibit sharp and hysteretic switching between superconducting and normal states with 
well-defined critical and retrapping currents. The two-dimensional (2D) critical current density 
e/2D,c is estimated to be as high as 1.8 A/m at 1.8 K. The temperature dependence of J2d,c indicates 
that the surface atomic steps play the role of strongly coupled Josephson junctions. 



PACS numbers: 68.35.B-,73.20.At,74.25.F- 

The state-of-the-art nanotechnology has enabled fabri- 
cation of ultrathin superconductors of high crystallinity 
and with atomically controlled thicknesses and interfaces 
1-|8|. This has opened ways to tune superconductivity 
l|, [sl, Isl and to investigate the thinnest crystalline lay- 
ers for its emergence [2|, |6|, |7| . The scope of the research 
on superconducting films has been substantially widened 
compared to the conventional studies where the samples 
were limited to granular and amorphous films [ol, llO| . 

Notably, superconductivity was found to exist for sil- 
icon surface reconstructions with metal adatoms [11|, 
which are the ultimate forms of thin epitaxial films. This 
is of primary importance because a variety of metal- 
adsorbed semiconductor surfaces [12| should now be re- 
garded as candidates for new superconducting materials. 
The above finding is, however, based on spectroscopic ev- 
idence of superconducting energy gaps observed by scan- 
ning tunneling microscopy (STM). Surface supercurrents 
could be very fragile, but such information was not avail- 
able. The local superconductivity shown there does not 
even guarantee the presence of macroscopic supercur- 
rent because a surface is inevitably severed by numerous 
atomic steps. The presence of macroscopic supercurrent 
can be best shown by electron transport measurements. 
Such a measurement has, however, been hampered by 
the requirements of an ultrahigh vacuum (UHV) and 
low-temperature environment and of stable contacts to 
surface atomic layers. The need to separate the surface 
current flow from the bulk contribution is also crucial 

mill- 

In this Letter, we perform direct and macroscopic elec- 
tron transport measurements on a silicon surface recon- 
struction with In adatoms (Si(lll)-( V7x V3)-In) [Tl',^- 
17| in UHV at low temperatures. The superconducting 
transition is evidenced by observations of the zero resis- 
tance state and of I — V characteristics exhibiting sharp 
and hysteretic switching below 2.8 K. This macroscopic 
supercurrent also shows a significant robustness; the two- 
dimensional (2D) critical current density J2d,c is esti- 



mated to be as high as 1.8 A/m at 1.8 K. The observed 
temperature dependence of J2d,c indicates that the sur- 
face atomic steps serve as strongly coupled Josephson 
junctions. 

Figure 1(a) shows an atomic structural model of 
Si(lll)-(V7 X Vsyin surface (referred to as {V7 x V^)- 
In) proposed by Kraft et al. [l5|, where one In atom 
corresponds to one unit cell of the Si(lll)-lxl surface. 
This surface hosts two-dimensional nearly free electrons 
with a high density [16|, exhibiting metallic transport 
properties down to low temperatures [17|. The central 
parallelogram indicates the \/l x \/3 unit cell. Note that 
the precise positions of In atoms have not been deter- 
mined and are assumed to be the hollow sites of the 
Si(lll)-lxl here. The surface was prepared by ther- 
mal evaporation of In onto a clean Si(lll) substrate fol- 
lowed by a brief annealing in UHV |15l4l7j , and was sub- 
sequently characterized by low energy electron diffraction 
(LEED) and STM. Sharp LEED spots (Fig. 1(b)) and 
high-resolution STM images (Fig. 1(c)) confirmed the 
presence of a well-ordered (>/7 x >/3)-In phase [l8|. For 
the large-scale morphology, STM observations show that 
the surface consists of flat terraces separated by atomic 
steps with a height of 0.31 nm (see Fig. l(d)(e)). The 
averaged terrace width was estimated to be 370 nm from 
repeated measurements. 

To perform four- terminal resistance measurements, the 
sample surface was patterned as depicted in Fig. 2(a). 
The central square with a size of 1 x 1 mm^ and the 
four outer squares connected at corners (purple areas) 
are made of the (\/7 x A/3)-In surface. The latter serves 
as current and voltage terminals for the former, thus en- 
abling van der Pauw's measurements [19|,l20|. The rest of 
the sample surface (gray areas) consists of bare Si, which 
is prepared by Ar+ sputtering throu gh a shadow mask 
(see Supplemental Materials (SM) 1) [2l|. Following the 
sample preparation, four Au-coated spring probes were 
brought into contact with the current/voltage terminal 
patterns in a UHV-compatible cryostat. Four-terminal 
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FIG. 1: (Color) (a) Atomic structural model of the (a/7x a/3)- 
In surface proposed by Kraft et al. ^IS]. The red and blue 
spheres represent In and surface Si atoms, respectively. The 
dashed lines indicate the {V7 x v^) unit cell. The In atoms 
are assumed to be on the hollow sites of the Si(lll)-lxl for 
simplicity, (b) LEED pattern taken at beam energy E — 86 
eV. The dashed lines indicate the reciprocal unit cell, (c) 
High-resolution STM image taken at a sample voltage Vs = 
—0.015 V and a tunneling current It = 270 pA. (d) Large- 
scale STM image taken at 14 = 2 V and It = 160 pA. (e) 
Height profile taken along the dashed line in (d). 
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FIG. 2: (Color) (a) Drawing of the sample patterning. The 
central and the four outer squares (purple areas) are made of 
(a/7 X a/3) -In, and the surrounding regions (gray areas) con- 
sist of sputtered Si surfaces, (b) Calculated current density 
distribution on the sample. The bright green (dark purple) 
represents a high (low) current density. The red dotted line 
indicates the flow of a bias current, (c) Temperature depen- 
dence of zero bias resistances. i?o,i (red lines) and i?o,ii (blue 
lines) are zero bias resistances measured using the probe con- 
figurations I and II, respectively, with DC bias currents of 1 
/iA. Left inset: i?o,i and i?o,ii for a larger temperature range. 
Right insets: schematic drawings of the probe configurations 
I and II. 



zero bias resistance Rq and I — V characteristics of the 
(a/7 X A/3)-In surface were measured from 20 to 1.8 K. 
Here Rq is simply defined as the ratio of the measured 
voltage to the bias current. The sample was carefully 
shielded from room temperature radiation and Si-diode 
thermometers were thermally anchored to prevent an er- 
roneous reading. To access the electron conduction only 
through the (a/7 x a/3) -In surface at low temperatures, 
Si(lll) substrates without intentional doping (resistiv- 
ity p > 1000 l^cm) were used [14]. The whole proce- 
dure including the pattern fabrication and the transport 
measurement was performed without breaking the UHV 
environment [IJ, |22| . 

Figure 2(c) exemplifies the temperature dependence of 
zero bias resistance. i?o,i and i^o,ii were measured with 
two complimentary configurations I, II by rotating the 
pairs of current/voltage probes by 90° (see the right in- 
sets of Fig. 2(c)). DC bias currents of 1 jaA were supplied 
and the offset voltages due to the thermoelectric effect 
were removed by inverting the bias polarity. Followed by 
gradual decreases with decreasing temperature T, i?o,i 
and i?o,ii dropped to nearly zero simultaneously below 
T = 2.8 K, indicating a superconducting transition. The 
residual resistances for T < 2.6 K are negligibly small 
compared to the noise level of 0.2 Q. The fact that both 
i?o,i and i?o,ii became zero excludes a possibility of fail- 



ure to detect a voltage drop due to an extremely high 
transport anisotropy. Although the transition is sharp 
just above the onset of the zero resistance state (ZRS), 
it exhibits a precursor below ~ 4 K which is absent for 
pure bulk superconductors. This is attributed to the su- 
perconducting fluctuation effects inherent to 2D super- 
conductors [23]. The transition temperature Tc = 2.8 K 
determined from the onset of the ZRS is slightly lower 
than the previously reported Tc = 3.14 K, which was 
determined from the opening of superconducting energy 
gap [HI. 

Assuming a homogeneous and isotropic resistivity for 
simpHcity, we calculated the current distribution within 
the present sample using the finite element method (Fig. 
2(b)). This allows us to determine the ratio P2b/Ro to 
be 4.54, where /)2D is the sheet resistance (2D resistivity) 
and Ro the four-terminal resistance defined as above. If 
(i?o,i + ^o,ii)/2 is identified with Rq^ p2D ^ 410 ft at 
5 K is obtained. This value is sufficiently smaller than 
the critical sheet resistance h/Ae^{= 6.45 x 10^ Q) of the 
superconductor-insulator transition, which is the crite- 
rion for the emergence of global superconducting coher- 
ence [Sl, llO[. Let us mention that the surface has ac- 
tually some anisotropy depending on samples and loca- 
tions. The origin of the anisotropy can be attributed to 



the local directions of surface steps, which should be the 
dominant electron scatterers at low temperatures. This 
may also cause the different temperature dependences of 
^0,1 5 ^0,11 displayed in Fig. 2(c). 

A further evidence for the presence of supercurrent was 
obtained by measuring I — V characteristics. The main 
panel of Fig. 3(a) shows a series of / — V characteristics 
taken by sweeping DC bias currents at different temper- 
atures from 3.11 K to 1.77 K. The bias current was swept 
toward the increasing direction at a rate of 87 /iA/s and 
the probe configuration I was adopted. The results were 
nearly independent of the sweeping rate between 23 and 
350 /iA/s. Below 2.8 K, the sample first switched from 
the normal state to the ZRS (where dV/dl = 0) at the 
retrapping current Ir and then from ZRS to the normal 
state at the critical current 1^ The switching behaviors 
became remarkably pronounced as the temperature was 
lowered down to 1.77 K. The ZRS can be safely assigned 
to the superconducting state because it is destroyed by 
an excessive current. By inverting the sweeping direction, 
the hysteresis oi I — V characteristics was confirmed as 
shown in the inset. The origin of the hysteretic switching 
can be a Joule heating effect ^]. Figure 3(b) summa- 
rizes the temperature dependence of critical current Ic 
(green squares) and retrapping current Ir (pink squares). 
The data taken with configurations I and II are shown 
by closed and open squares, respectively. Ic and Ir were 
almost identical for the two configurations. Following the 
onset around 2.8 K, both Ic and Ir steadily increase as 
temperature is lowered, reaching 520 /iA and 230 /iA at 
1.8 K, respectively. 

The 2D critical and retrapping current densities J2D,ci 
J2D,r are determined as described below [21]. First we 
note that the self-screening of magnetic field produced by 
supercurrent is very weak because the superconducting 
layer is atomically thin in the present system (SM 2). In 
addition, if the magnitude of the superconducting order 
parameter is constant, the distribution of supercurrent is 
identical to that of normal current when the same bound- 
ary condition is imposed (SM 3). Therefore, the above 
calculation on the normal current density is also valid 
here. Figure 2(b) shows that the current density is the 
highest at the constrictions between the current probes 
and the central area, which means that Ic is determined 
in the constrictions (SM 4). Since the current density 
is nearly constant in the middle of the constriction, the 
measured Ic and Ir can be converted to J2d,c and J2D,r 
by dividing them by its width Wc = 0.283 mm (see the 
right axis of Fig. 3(b)). J2d,c = 1-8 A/m at 1.8 K is 
remarkably high considering that the conducting layer is 
single- atom thick. If the thickness of (V^ x ^/3)-ln is 
assumed to be double the covalent radius of In (= 0.30 
nm), this corresponds to a 3D critical current density 
J3D,c = 6.1 X lO^A/m^. 

The mechanism of determining the critical current is 
discussed as follows. At low temperatures, the atomic 



> 

E 

CD 

B 
O 

> 




3.11 K 

2.73 K Iq_| 

2.61 K 

2.51 K 

2.20 K 

2.01 K 

1.77 K 



-600 -400 -200 200 400 600 
Bias Current (|liA) 



< 




2.0 



- 1.5 



- 1.0 S 



0.5 



2.0 2.4 2.8 

Temperature (K) 



3.2 



0.0 



FIG. 3: (Color) (a) Temperature dependence of I — V charac- 
teristics obtained with configuration I. The inset shows I — V 
characteristics taken by inverting sweeping directions at 1.80 
K. The arrows indicate the sweeping directions, (b) Temper- 
ature dependences of critical current Ic (green squares) and 
retrapping current Ir (pink squares). The data taken with 
configurations I and II are shown by closed and open squares, 
respectively. The blue solid and red dotted lines show theo- 
retical fits. For details, see the text. The inset sketches the 
relation between current flow and atomic steps. 



steps are considered to be the dominant source of elas- 
tic electron scattering and resistance on the surface. In 
the superconducting state, therefore, they can serve as 
Josephson junctions. In this case, the temperature de- 
pendence of J2D,c is given by the following equation [ 25| : 



J2dAT) 



ttAJT) 



tanh{A{T)/2kBT), 



(1) 



where A(T) is the superconducting energy gap at T, 
Pstep the (T-independent) normal resistance of the atomic 
step for unit length, ks the Boltzmann constant. Since 
A(T) follows the BCS theory ^ and T^ = 2.8 K is al- 
ready known, the experimentally obtained J2D,c{T) can 
be fitted with Eq. ([!]) using J2d,c(0) as the only param- 
eter. The fitting can reproduce the experiment semi- 
quantitatively, which gives J2d,c(0) = 2.7 A/m (the 
blue solid line in Fig. 3(b)). Another requirement of 



A(0) = 0.57 meV [ll| allows us to determine pstep to be 
3.3 X 10~^ Qm. If the critical current is determined by 
Cooper pair breaking induced by an excessive supercur- 
rent, it follows the equation J2d,c(T) cx (1 - [T/Tcff^'^ 
\Z&\. Fitting with the equation gives a poor result (red 
dotted line in Fig. 3(c)), excluding this mechanism. This 
suggests that the supercurrent density is not high enough 
to substantially suppress the magnitude of order param- 
eter [27], in consistent with the earlier assumption on its 
spatial uniformity. We note that the surface steps should 
be regarded as strongly coupled junctions because of the 
high critical currents, although the term Josephson junc- 
tion is conventionally used for a weak coupling. 

The normal resistance of atomic steps pstep obtained 
above can be compared to that from the normal sheet re- 
sistance P2D • Suppose a local current flows across atomic 
steps at an angle of and the steps are separated by an 
average distance d (see the inset of Fig. 3(b)). If we 
simply assume that is randomly distributed between 
-7r/2 and 7r/2 over the sample surface, the average sheet 
resistance p2D can be calculated as 



P2D 



7r/2 
-7r/2 



Pstep I sin l9 1 _ 2pstep 



d 



nd 



(2) 



Insertion of experimentally obtained values d = 370 nm 
and P2D = 410 ft into Eq. ([2j) gives pstep = 2.4x 10~^ ftm. 
This is in satisfactory agreement with the value pstep = 
3.3 X 10~^ ftm determined earlier, supporting the Joseph- 
son junction model of atomic steps. We note that, in the 
analysis of the Josephson junction, the angle 6 between 
the current flow and steps was not taken into account. 
This is because the critical current is determined by indi- 
vidual atomic steps (presumably by one with the highest 
Pstep) and the step separation d is not relevant. On the 
contrary, in the case of sheet resistance, a smaller in- 
creases the effective step separation as d/sin^, resulting 
in a lower p2D- We also note that the obtained step 
resistances are comparable to the previously reported 
Pstep ^ 2 X 10-"^ ^m for Si(lll)-(V3 x V^)-Ag [28j, 
although a different surface reconstruction was studied 
there. 

In conclusion, we have demonstrated that macroscopic 
and robust supercurrents can run on a {^/7 x ^/3)-In sur- 
face despite the presence of atomic steps. It was indicated 
that the surface atomic steps serve as strongly coupled 
Josephson junctions. The present study makes various 
surface reconstructions of silicon and related semiconduc- 
tors candidates for practical superconducting materials. 
We envision that atomic-scale design and tuning of su- 
perconductivity will be feasible for such surface systems 
based on the current nanotechnology. 
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